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Levy ratchets are minimal models of fluctuation-driven transport in the presence of Levy noise 
and periodic external potentials with broken spatial symmetry. Levy noise can drive the system 
out-of thermodynamics equilibrium and previous works have shown that in this case, a net ratchet 
current can appear even in the absence of time dependent perturbations, external tilting forces, or 
a bias in the noise. The majority of studies on the interaction of Levy noise with external potentials 
have assumed a-stable Levy statistics in the Langevin description, which in the continuum limit 
corresponds to the fractional Fokker-Planck equation. However, the divergence of the low order 
moments is a potential drawback of a-stable distributions because, in applications, the moments 
represent physical quantities. For example, for a < 1, the current (first moment), J, in a-stable Levy 
ratchets is unbounded. To overcome this limitation, in this paper we study ratchet transport using 
truncated Levy distributions which in the continuum limit correspond to the spatially tempered 
fractional Fokker-Planck equation. The main object of study is the dependence of the ratchet 
current on the level of tempering, A. For finite tempering, A ^ 0, the statistics ultimately converges 
(although very slowly) to Gaussian diffusion in the absence of a potential. However, it is shown here 
that in the presence of a ratchet potential a finite non-equilibrium current persists asymptotically 
for any finite value of A. The current is observed to converge exponentially in time to the steady 
state value. The steady state current exhibits algebraically decay with A, as J ~ A , for a > 1.75. 
However, for a < 1.5, the steady state current decays exponentially with A, as J ~ e~^ A . In the 
presence of a bias in the Levy noise, it is shown that the tempering can lead to a current reversal. 
A detailed numerical study is presented on the dependence of the current on A and the physical 
parameters of the system. 



I. INTRODUCTION 



The problem of fluctuation driven transport has a long 
and interesting history, with early discussions dating 
back to work by Smoluchowsky in 1912. For a review, 
see Ref. pQ and references therein. Of particular interest 
are periodic potentials with broken spatial symmetry, 
known as ratchet potentials, see for example Fig. [T] In 
this case, particles in the potential wells experience an 
asymmetric force and one might naively expect that 
in the presence of thermal equilibrium fluctuations 
particles will, on the average, drift down-hill in the 
direction of the weaker force, giving rise to a net current. 
However, it has long been recognized that this naive 
picture is incorrect since net work cannot be extracted 
out of equilibrium thermal fluctuations. However, in the 
presence of non-equilibrium perturbations the situation 
is quite different and net currents can be generated. 
For example, in the case of pulsating ratchets a current 
can be generated by driving the system out of thermal 
equilibrium by adding a time variation to the amplitude 
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of the potential. Considerable amount of work has thus 
been devoted to construct "minimal" models of ratchet 
transport which involve various far from equilibrium 
mechanisms. However, most of these studies have 
focused on models in which the fluctuations satisfy 
Gaussian statistics. 

The early work reported in Refs. [3] originally pro- 
posed a minimal model for Levy ratchets consisting of 
a time-independent ratchet potential in the presence of 
uncorrelated Levy noise. It was shown that, even in the 
absence of an external tilting force or time dependence 
in the potential, or a bias in the noise, the Levy flights 
drive the system out of thermodynamic equilibrium 
and generate an up-hill current (i.e., a current in the 
direction of the steeper side of the asymmetric poten- 
tial). Following this, Ref. [4] studied the Levy ratchet 
problem and proposed robust probability measures of 
directionality of transport. In Ref. |5j the influence 
of periodically modulated Levy noise asymmetry was 
studied. Beyond the study of ratchet transport, the 
study of the interaction of Levy noise with external po- 
tentials has been a topic of considerable recent interest. 
Some examples include the study harmonic [6l [7] and 
non-harmonic |8HT0] potentials, and the study of the 
Kramers' problem in Refs. [TTMl4| . 
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Levy statistics has been used to model a wide range of 
problems involving probability density functions (pdfs) 
with slowly decaying tails. However, an often overlooked 
issue is the differences that typically exist between the 
strict mathematical properties of Levy distributions and 
the pdfs found in practical problems involving experi- 
ments or numerical simulations. For example, whereas 
the second and higher moments of a-stable Levy distri- 
bution with 1 < a < 2 are formally infinite, the moments 
of pdfs of practical interest are finite. One of the reasons 
why pdfs of practical interest have finite moments is that 
they typically have finite support because the statistical 
samples upon which they are constructed are finite. For 
example, pdfs with algebraic decaying tails have been ob- 
served in the description of particle displacements in tur- 
bulent transport [15] . These numerically determined pdfs 
have finite moments because their support is limited by 
the largest possible displacement which is finite. Another 
key issue is that in applications, the moments represent 
physical quantities that cannot be infinite. This issue is 
particularly delicate in the study of Levy flights in ex- 
ternal potentials where a-stable Levy noise is added to 
the Langevin stochastic equation for the particle veloc- 
ity. In this case, the divergence of the second moment can 
lead in principle to a divergence in the kinetic energy. In 
practical terms this problem might not be disastrous be- 
cause any numerical realization of the numerical model 
will have finite moments. Nevertheless this shortcom- 
ing of the Levy a-stable Langevin model, which is also 
present in the corresponding a-stable fractional Fokker- 
Planck model, has to be resolved. In this paper we pro- 
pose to use the truncated fractional diffusion operators 
introduced in Ref. |16j to address this problem. 

Truncated Levy distributions distributions were orig- 
inally introduced in Ref. |17l 118] as a simple prescrip- 
tion to guarantee the finiteness of the second moment. 
Reference [E5] considered a general class of multivari- 
ate tempered stable Levy processes and established their 
parametrization and probabilistic representations. Gen- 
eral Levy processes, and in particular tempered stable 
processes, were incorporated in the continuous time ran- 
dom walk model and in macroscopic transport models 
in Ref. [16]. This reference introduced the notion of 
truncated fractional derivatives in the construction of 
non-diffusive transport models driven by truncated Levy 
flights. The role of tempered Levy distribution in the 
super-diffusive propagation of fronts in reaction-diffusion 
systems has been recently studied in Ref. [2"0] , 

The goal of this paper is to study the role of truncation 
on fluctuation driven transport with Levy noise. Our ap- 
proach is based on the numerical solution of the spatially 
tempered Fractional Fokker Planck equation, obtained 
by replacing the diffusion operator by a truncated frac- 
tional diffusion operator. We use two complementary 
numerical methods. One is based on the finite differ- 
ence, Grunwald-Letnikov discretization of the truncated 
fractional derivatives that allows the computation of the 



space-time evolution of the pdf in a finite size domain. 
The second method is a Fourier based spectral method 
that allows the computation of the reduced pdf in a pe- 
riodic domain as well as the time dependent and asymp- 
totic steady state current. The main object of study is 
the dependence of the steady state current on the level of 
truncation, A, the stability index a, and the noise asym- 
metry 9 as well as the asymmetry of the potential. 

The rest of the paper is organized as follows. In the 
next section we discuss the setting of the Levy ratchet 
problem and introduce the spatially tempered Fractional 
Fokker Planck equation. Section III presents the numer- 
ical results obtained from the solution of the time depen- 
dent and the steady state tempered Fractional Fokker 
Planck equation. Section IV contains the conclusions. 

II. SPATIALLY TEMPERED FRACTIONAL 
FOKKER-PLANCK EQUATION 

The most natural formulation of the problem of 
fluctuation-driven transport in the presence of an exter- 
nal potential is based on the Langevin approach. The 
starting point is an ensemble of particles with coordinates 
x(t) which in the over-damped case evolve according to 
the stochastic differential equation 

dx = - (d x V) dt + dr] (1) 

where for simplicity, and without loss of generality, we 
have assumed a one-dimensional domain. V(x) denotes 
the external potential, and drj is a stochastic processes 
modeling the fluctuations. In the case of uncorrelated, 
Gaussian fluctuations, drj — y/2~xdW, where dW is a 
stochastic Wiener process, with (dW) — and (dW 2 ) — 
dt. As it is well-known, if the individual particles evolve 
according to the Langevin equation (JlJ with uncorre- 
lated Gaussian noise, then the probability density func- 
tion (pdf) satisfies the Fokker-Plank equation 

d t P = d x [Pd x V]+xd 2 x P. (2) 

Of particular relevance to the present paper is the 
study of periodic potentials with broken spatial sym- 
metry. These, potentials, also called "ratchets", satisfy 
V{x) = V(x + L) where L is the length of the period but 
they lack reflection symmetry. A paradigmatic example 
is V = Vb[sin(27rx/L) + 0.25sin(47ra;/L)]. For an easier 
control of the degree of spatial symmetry, following [3J, 
here we use 

y _ y ( 1 — cos [nx/ai] if < x < a\ , , 

1 1 + cos [tt(x — ai)/a 2 ] if ai < x < L , ^ ' 

where Vq is the amplitude, L — a\ + a 2 is the period, 
and A — (a± — a 2 )/L is the asymmetry parameter. In 
all the calculations presented here, Vq = L = 1. In the 
study of ratchets it is customary to add an external con- 
stant "tilting" force F to the potential, and to consider 
the effective potential V e ff = V — Fx. However, unless 
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FIG. 1: Ratchet potential in Eq. Q for A = -0.5 and V = 1. 

mentioned otherwise, the calculations presented here will 
for the most part assume F = 0. Figure [T] shows a plot of 
the ratchet potential in Eq. Q with a% = 1/4, a% = 3/4, 
Vb = 1, and L = 1, that correspond to A = —0.5. 

The conditions for the existence of ratchet currents 
in the presence of diffusive transport have been exten- 
sively studied and are well-understood, see for example 
Ref. [T] and references therein. However this is not the 
case for fluctuation-driven transport in the presence of 
non-diffusive transport. To study the role non-diffusive 
transport in fluctuation-driven ratchet transport, Ref. [3J 
considered the Langevin model but with dr\ correspond- 
ing to a Levy process. In this case, the continuum limit 
of the Langevin equation leads to the a-stable fractional 
Fokker-Planck equation 



where 



dtP = d x [Pd x V]+xd%P. 



d«p = i_ 00 D™p + r x D«p. 



(4) 



(5) 



The operators on the right hand side of Eq. ^ are the left 
and right Riemann-Liouville fractional derivatives which 
are integro-differential operators defined in Fourier space 
according to [2"HI23] 



T [-coD-P] - (-»fc) Q P , T [ X D^P] = {ik) a P , 



where F[f] (k) = f — J f(x)e lkx dx. The factors 



(I- 



2cos(a7r/2) 



(1 + 0) 
2cos(mr/2) ' 



(6) 



(7) 



determine the relative weight of the left and the right 
fractional derivatives in terms of the parameter 6. The 
order of the fractional operators a, the asymmetry pa- 
rameter 6, and the diffusivity \ correspond to the index, 
the skewness, and the scale factor of the characteristic 
function of the Levy noise in the Langevin formulation. 

From the statistical mechanics point of view, the dif- 
ference between the transport operators in Eqs. ([2| and 



Q rests on the different assumptions of the underlying 
stochastic process. As it is well known, the standard 
diffusion operator in Eq. |2]) assumes a Brownian ran- 
dom walk, whereas, in the context of the continuous time 
random walk (see for example Ref. and references 
therein) the fractional derivative operators in Eq. Q 
assume an a-stable Levy process describing the jumps 
of particles. One might naively expect that these two 
cases encompass all the fundamentally different trans- 
port operators in the absence of memory and correla- 
tions. This expectation is based on the fact that accord- 
ing to the generalized central limit theorem the sum of 
individual random particle displacements asymptotically 
converge to either a Gaussian (described by the stan- 
dard diffusion equation) or an a-stable Levy distribution 
(described by the a-stable fractional diffusion equation). 
However, an issue of significant physical relevance is that 
the convergence to the Gaussian or a-stable behaviors 
could be extremely slow in problems of practical inter- 
est. For example, in the case of truncated Levy flights, it 
is known [17l[24l|25] that the crossover time scale, r c , of 
the convergence to Gaussian behavior scales algebraically 
t c ~ \~ a where A is the level of truncation. Such slow 
algebraic convergence indicates that for the intermedi- 
ate time scales of practical interest, these type of pro- 
cesses cannot accurately be described as Gaussian or a- 
stable. Motivated by this, Ref. [T5] considered the contin- 
uous time random walk for jump processes corresponding 
to general Levy processes in the Levy-Khintchine repre- 
sentation, and obtained integro-differential operators de- 
scribing transport in the continuous, fluid limit. These 
operators contain as special cases the regular (Laplacian) 
and the a-stable fractional diffusion operators (Riemann- 
Liouville fractional derivatives), and, most importantly, 
allow the incorporation of a richer class of stochastic pro- 
cesses of physical relevance. In particular, for the case of 
truncated Levy processes, Ref. [16] proposed the trun- 
cated fractional diffusion transport operator 



d"' x P = V"' X P + vd x P - vP , 



(8) 



where T> X ' X is the A-truncated fractional derivative oper- 
ator of order a, defined as 



V a - X = le- Xl 



,££ e Xx + re Xx X D^ e~ Xx , (9) 



where -oqD" and X D^ are the Riemann-Liouville deriva- 
tives in Eq. Q. The effective advection velocity is de- 
fined as 



and 



< a < 1 

|cos(a7r/2)| 



A" 



cos (a7r/2) 



(10) 



(11) 



As expected, for A = 0, the transport operator in Eq. ^ 
reduces to the a-stable fractional diffusion operator in 
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Eq. ([5]). Based on the previous discussion, we pro- 
pose the following spatially-tempered Fractional Fokker- 
Planck equation 
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d t P = d x [Pd x V]+ X d^ x P : 



(12) 



to study the role of truncation effects on Levy flights in 
the presence of external potentials in the intermediate 
asymptotic regime. In particular, the study presented 
here of the role of truncation in fluctuation-driven trans- 
port in the presence of Levy flights is based on the numer- 
ical integration of Eq. ( 12 ). Taking the Fourier transform 
of Eq. (pi), we get 



d t P = A(k)P - ikT[Pd x V] (k) , 



(13) 



where A(k) is the logarithm of the characteristic expo- 
nent of the tempered Levy process 



A = 



_ -x 



2cos(a7r/2) 



(1 + 0)(A + ik) a + (1 - 0)(A - ik) a - 2A Q , 

(1 + 0)(A + ik) a + (1 - 0)(A - ik) a - 2\ a - 2ika6X a - 

(14) 

for < a < 1 and 1 < a < 2 respectively. 



III. NUMERICAL RESULTS: 
CURRENT 



RATCHET 



The spatio-temporal evolution of the pdf in the finite 
domain (0, 1) was obtained by solving numerically the 
spatially tempered fractional Fokker-Planck equation us- 
ing a finite difference method. Following Ref. |26j . we use 
regularized (in the Caputo sense) fractional derivatives in 
space. However, here we factorize the regularized opera- 
tors as %D a £ =§ D x {2 ~ a) d 2 x , and C X D? =% D~ {2 ~ a) d 2 x , and 
discretize the fractional integral operators, gD^ ^ ^ and 



id; 



(2-0) 



using the Grunwald-Letnikov representation. 



The second order derivative, d x , was discretized using 
a central difference method. Following Ref. [3], we dis- 
cretized the time evolution using an operation splitting 
method that separates the fractional derivative operators 
from the advection and potential terms. The fractional 
derivative operators (half) time step was done using a 
weighted average Crank-Nicholson method. 

Figure [2] shows snapshots in time of the numerically 
calculated pdf's for different levels of truncation. The 
initial condition corresponds to a distribution of parti- 
cles localized in the potential well at the middle of the 
computational domain 



P(x,0) = 



'2tt 



exp 



(X-1/2Y 



(15) 



with a = 0.12. Compared to the A = (a-stable) case, 
it is observed that the tempering reduces the "leakage" 
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FIG. 2: Time dependent solutions of the truncated factional 



Fokker-Planck equation ( 12 1 for A = and A = 3, with a = 



1.25, 6 = 0, Vb = 1, and \ = 0.5. The solid line at the bottom 
is the potential in linear-linear scale. The pdfs for A = and 
A = 3 are sown in a log-linear scale. 



of the pdf out of the potential well. However, for both 
A = and A = 3 the profile peaks are higher on the right 
than on the left, indicating the presence of a net current. 
The existence of the current for A = is consistent with 
the results previously reported in Ref. [3J. The under- 
standing and characterization of the current for A 7^ is 
one the objectives of this section. 

Given a solution, P, of the tempered fractional Fokker- 
Planck equation, the current is defined as the rate of 
change of the first moment, 



J = 




k=0 



From Eqs. (16) and (13), it follows that 



J=-T[Pd x V] k=Q 



dA 
~d~k: 



k=0 



In the a-stable case, 



A = 



dA 
dk 



k=0 



00 





< a < 1 
1 < a < 2. 



(16) 



(17) 



(18) 



According to this, without tempering, the current is not 
well-defined for a < 1 due to the divergence of the first 
moment of Levy distributions. However, the introduction 
of tempering leads to a well-defined current because in 
this case 



A ^ 



dA\ 
dk ) k=0 



|cos(q7 



/2)l 





< a < 1, 
1 < a < 2 



(19) 



The asymptotic steady state current can be obtained 
numerically by advancing the solution for P to long times 
as done in Ref. [3] . However, here we follow a more direct 
approach based on the computation of the steady state 
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solution in Fourier space, which can be formulated as 
an eigenvalue problem as shown below. For a general 
periodic potential, V{x + L) = V(x), 



d x V= £ 

n— — oo 

where k n — 2-kti/L, we have 

oo 

F[Pd x V]= c n P(k + k n ,t) 



(20) 



(21) 



Evaluating this expression at k — and substituting into 
Eq. (17) we get the following expression for the current 



J 



V c P 



dA 
dk 



(22) 



k=0 



where P n — P (k = k n ) is the Fourier transform of the 
steady state solution. To determine the time ev olut ion 
of the n-th Fourier mode, P n (t), we evaluate Eq. (131 at 
k = k n and get the following set of coupled differential 
equations 



OtPn 



A„P„ 



+3 ' 



(23) 



j = -oo 



where A n = A(fc = k n ). The steady state asymptotic cur- 
rent, J(oo), is obtained by setting d t P n — in Eq. (23) 
and solving the set of algebraic equations 



j r n+j 



0. 



(24) 



j = -oo 



For the numerical solution of Eqs. (23) and (24), we 
truncated the infinite series at a large value of n — N 
and solved the system with the normalization constrain 
Pq = 1 (as required by the conservation of probability) . 
In the calculations reported here N ranged from 500 to 
2000. The unique solution {P n } was then used in Eq. ( 22 ) 
to evaluate the current. The corresponding reconstructed 
pdf in real space gives the reduced pdf, P r , 



N 



(25) 



-A' 



which is the steady state solution of the tempered frac- 
tional Fokker-Planck equation in a periodic domain. This 
solution can also be obtained by "folding" the unbounded 
domain solution in the (0, L) domain pQ 



P r = E P ( x + nL > l ) 



(26) 



Figure [3] shows the reduced steady state pdf in the 
a-stable case, i.e., A = 0, for different values of a. 
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FIG. 3: Steady state solution of the reduced pdf in the pe- 
riodic domain for different values of a, and A = 0, 9 — 0, 
V = 1,A = -0.274, and X = 0.5. 




FIG. 4: Steady state solution of the reduced pdf in the pe- 
riodic domain for different values of A, and a — 1.25, 6 — 0, 
Vo = 1, A = -0.274 and x = 0.5. The solid line is the ratchet 
potential. 



In all the numerical results, the current has been non- 
dimensionalized using the scale L/r, where r = L a /x 
is the fractional diffusion time scale. As expected, for 
a = 2, P r corresponds to the Boltzmann distribution 
P r {x) = (1/Z) exp (— V{x)/x), where Z is the normaliza- 
tion constant. As the value of a decreases, the pdf sig- 
nificantly departs from the Boltzmann distribution and 
develops the asymmetry responsible for the finite net cur- 
rent observed for a < 2. The dependence of the reduced 
pdf on A for a fixed value of a is shown in Fig. [4] The leak- 
age of the pdf outside the potential wells is significantly 
reduced for large values of A. As we will discuss later, in 
general the current decreases with increasing A, but the 
pdf does not approach the diffusive Boltzmann distribu- 
tion. Figure [5] shows the dependence of the current on 
the value of a for different levels of tempering. In agree- 
ment with the results reported in Ref. 3J , in the absence 
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FIG. 5: Ratchet current, J, as function of a, for different 
levels of truncation, A, with 9 = 0, A = —0.6, Vo = 1, and 
X = 0.5. 



of tempering, A = 0, a net current is observed for a < 2 
with a maximum around a » 1.4. Most importantly, as 
the hgure shows, a finite (albeit smaller) current remains 
for A ^ 0. That is, although for A ^ 0, the statistics 
of the tempered fractional operator in Eq. ^ ultimately 
converges to Gaussian behavior [T^] , in the presence of a 
ratchet potential a net non-equilibrium current persists 
in the steady state asymptotic regime. As expected, re- 
gardless of the value of A, in the limit a = 2, the noise 
becomes Gaussian and the current vanishes. It is ob- 
served that as the truncation increases, the value of a 
for which the maximum current is attained shifts to the 
right, and for large values of A the strength of the current 
exhibits only a weak dependence on a. To study the rate 
of convergence to the steady state asymptotic current, 
we show in Fig. [6] S J(t) = J(t) — J(oo) as function of 
time where J(i) is the time dependent transient current 
and J(oo) is the steady state asymptotic current. An 
exponential convergence of the form 



5 J 



,-vt 



(27) 



is observed, and as Table I shows, the decay rate, 77, 
exhibits a very weak dependence on the value of a and 
A. 



TABLE I: Dependence of exponential decay rate, 77, of tran- 
sient current in Eq. ( |27| l as function of a and A for 9 — 0, 
A = -0.274, Vo = 1, and x = 0.5. Case a = 1.5, A = 10 is 
shown in Fig. [6] 
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1.25 
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64 


58 


58 


60 


57 


58 


52 


51 


53 



As the previous calculations show, as A increases the 
current is reduced, and a problem of significant interest is 



FIG. 6: Exponential convergence of ratchet current as func- 
tion of time. The plot shows 8 J = J (pa) — J(t) where J(oo) 
is the steady state asymptotic current (obtained from the so- 
lution of Eq. (24l) and J(i) is the time dependent transient 
current (obtained from the solution of Eq. (23 1. The param- 
eter values are a = 1.5, X = 10, 9 = 0, A = -0.274, Vo = 1, 
and x = 0.5. The dashed line is an exponential fit of the form 
8 J ~ e~ vt with 77 ~ 57. Table I shows the exponential decay 
rate for different values of a and A. 



to find the rate of decay of the current in the asymptotic 
limit A — > 00. As Figs. 7]|9 show, depending on the value 
of a, two asymptotic regimes are observed. For a near the 
Gaussian value a = 2, the current exhibits an algebraic 
decay for large A of the form 



J ~ A~ ? . 



(28) 



where £ w 2.5 for a — 1.9 and C 1=3 6 for a = 1.75. This 
slow algebraic scaling breaks down as a is reduced from 
the Gaussian value. In particular, as Fig. [8] shows, for 
1 < a < 1.5 the current decays exponentially fast 



J 



-a 



(29) 



with £ « 0.40. As Fig. [9] verifies, the same scaling is 
observed for a < 1. 



Figure [10| shows the dependence of the current, J, on 
the potential asymmetry, A. As expected, in all cases the 
current vanishes when the potential is symmetric, A — 
0. Consistent with the results reported in [3], the curve 
A = 0, which corresponds to an a-stable Levy ratchet, 
shows the existence of a ratchet current in the direction 
of the steepest side of the potential. The main result of 
this study is that, even in the presence of truncation, a 
ratchet current is observed and its dependence on A is 
qualitatively similar to the a-stable case. To study the 
dependence of the current on an external tilting force, we 
computed the current with an effective potential, V e ff = 
V(x) — Fx where V(x) is the ratchet potential in Eq. (jjjj, 
and F is a constant. As shown on Fig. [TT] the A = 
case recovers the a-stable Levy ratchet results, and, as 
A increases, the strength of the current decreases. Note 
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FIG. 7: Decay of ratchet current, J, as function A for different 
values of a and 9 = 0, V = 1, and A = -0.274. The Log-Log 
scale shows evidence of algebraic decay for a > f .75. 



A 

FIG. 10: Ratchet current, J, as function of ratchet potential 
asymmetry, A, for different levels of truncation, A, a = 1.5, 
6 = 0, Vo = 1, and x = 0.5. 




FIG. 8: Decay of ratchet current, J, as function A for different 
values of a and 9 = 0, Vb = 1, and A = -0.274. The Log- 
Linear scale shows evidence of exponential decay for a < 1.5. 




FIG. 11: Ratchet current, J, as function of external force, F, 
for different levels of truncation, A, a — 1.5, 9 = 0, Vo = 1, 
4 = -0.274, and x = 0.5. 
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FIG. 9: Same as Fig. H] but for a < 1. 



that the value of the stopping force, F ps — 1, i.e. the 
force needed to cancel the ratchet current, seems to be 
independent of A. 

An interesting effect of the truncation is observed when 
the Levy noise is asymmetric, i.e. when the weighting fac- 
tors I, r of the left and the right fractional derivatives in 
Eq. @ are different. As shown in Fig. [12] in the absence 
of truncation, A = 0, the current can be reversed by bi- 
asing the asymmetry of the noise. In particular, for the 
case shown, increasing 9 to a large enough value stops the 
positive current. This is consistent with the fact that for 
6 > the skewness of the Levy noise is negative and the 
probability of having long jumps to the right is reduced. 
As the value of 9 is further increased the probability of 
jumps to the left increases and this eventually leads to a 
negative current. A similar phenomenology is observed 
when the truncation is present, except that (in addition 
to the usual overall decrease of the magnitude of the cur- 
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FIG. 13: Ratchet current, J, as function of truncation param- 
eter A, for different levels of fractional diffusion asymmetry 9, 
a = 1.5, Vo = 1, A = -0.274, and X = 0.5. For 6 = 0.275, 
the current is reversed for A « 1/3. 



rent) the critical value of 9 for current reversal increases. 
The fact that the critical 9 for current reversal depends 
on A implies that for a range of 9 values it is possible to 
reverse the current by changing A only. This interesting 
feature of truncated Levy ratchets is clearly illustrated in 
Fig. |l~3| where it is shown that for 9 = 0.275 the trunca- 
tion can in fact lead to a current reversal. This behavior 



is related to the fact that, as Eq. (10 1 shows, for 9 ^ 



the truncation gives rise to a finite advection term in the 
Fokker-Planck equation. To conclude the numerical sim- 
ulations of the steady state current, we show in Fig. [14] 
the dependence of the current on a = y^l a , where \ is 
the diffusivity. In the absence of truncation the current 
is maximum for a ~ 1. However, as the truncations in- 
crease the value of the diffusivity for which the maximum 
current is attained increases. 



IV. CONCLUSIONS 

The incorporation of tempering in models of anoma- 
lous transport due to Levy flights is motivated by the 
fact that in problems of practical interest, the moments 
of the pdfs represent physical quantities that cannot be 
unbounded. For example, although anomalously large 
displacements have been well documented in numerical 
simulations of turbulent transport and in fluid mechanics 
experiments, it is clear from the physical point of view 
that particle displacements cannot be arbitrarily large 
and that a cut-off or truncation has to be incorporated 
in the pdfs describing these processes. This was the intu- 
ition behind the original proposal of truncated Levy flight 
in Ref. [17] and the main motivation for the construc- 
tion of spatially tempered fractional diffusion operators 
in Ref. [TB]- In this paper we have studied the effects of 
spatial tempering on directed transport in Levy ratchets. 

Our study was based on the numerical solution of the 
time dependent and the steady state spatially tempered 
fractional Fokker-Planck equation. We used two com- 
plementary numerical methods. For the computation of 
the space-time evolution of the pdf in a finite size do- 
main, we used a finite difference method based on the 
Grunwald-Letnikov discretization of the truncated frac- 
tional derivatives. For the computation of the reduced 
pdf in a periodic domain as well as the time dependent 
and asymptotic steady state current, we used a Fourier 
based spectral method. The main object of study was 
the dependence of the steady state current on the level of 
truncation, A, the stability index a, and the noise asym- 
metry 9 as well as the asymmetry of the potential. 

The general conclusion is that in most parameter 
regimes the inclusion of truncation does not change the 
qualitative dependencies of the current on the different 
system parameters; instead, the effect of truncation is 
mainly to decrease the magnitude of current. However, 
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an interesting exception to this rule occurs when the Levy 
noise is allowed to have a bias. In this case, truncation 
can lead to a reversal of current direction. The persis- 
tence of ratchet currents for A ^ is interesting since the 
finite second moments of the truncated Levy statistics 
guarantee an eventual convergence to Gaussian statistics 
(in the absence of an external potential) and, and as it 
well-known, Gaussian fluctuations cannot create a net 
current in the absence of time dependent perturbations 
or an external tilting force. The current is observed to 
converge exponentially in time to the steady state value, 
and the decay rate exhibits only a weak dependence of 
a and A. The decay of the current for increasing values 
of A exhibits different scaling properties depending on 
the value of a. For large values of alpha, 1.75 < a, the 



current exhibits an algebraic decay and for small values, 
a < 1.5, the current exhibits an exponential decay. 
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